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Abstract 

We  show  that  a  strong  laser  acting  on  an  electron  in  a  quantum  double  well 
which  is  under  the  influence  of  a  strong  static  field,  will  cause  (for  special  values 
of  the  parameters)  emission  of  intense  low  frequency  radiation  and  can  either 
completely  delocalize  or  completely  localize  the  electron. 
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The  response  of  an  electron  in  a  double  well  (or  that  of  a  two  level  system 
representing  a  double  well)  to  a  strong  laser  has  been  the  subject  of  several  recent 
investigations10.  The  parameters  most  useful  for  describing  the  observations  are, 
e0  =  2£/hco  and  the  frequency: 

Qq  =  (2e/fco>)  J0(e0)  (1) 

Here  2e  is  the  energy  gap  between  the  lowest  two  levels,  co  is  the  laser  frequency, 
E0  is  the  electric  field  of  the  laser,  |i12  is  the  transition  dipole  between  these  levels 
and  J0  is  the  Bessel  function  of  order  zero11.  The  values  of  {o),E0}  for  which 
Qq  =  0  (2) 

are  called6"7  points  of  accidental  degeneracy  (AD  points).  According  to  (1)  they 
satisfy 

e0  =  2ni2E(/ho)  =  zi(0),  (3) 

where  Zj(0)  denotes  the  i-th  zero  of  J0.  If 

2e/hco  <1,  (4) 

the  behavior  of  the  system  when  {co,E0}  is  at  or  near  an  AD  point  is  very 
interesting.  As  is  well  known,  if  at  a  given  moment  the  electron  is  localized  in 
one  of  the  wells,  it  will  oscillate  between  the  wells  as  the  time  evolves.  A  laser 
with  the  parameters  chosen  at  an  AD  point  [i.e.  satisfying  Eq.(3)]  can  stop  this 
oscillation  and  keep  the  electron  in  one  of  the  wells1,2.  Moreover,  the  system  can 
emit  (or  absorb)  intense  even-harmonics6,7  [i.e.  the  induced  dipole  p(t)  has  intense 
Fourier  components  at  the  frequency  2nco,  where  n  is  an  integer].  The  presence  of 
this  emission  is  remarkable  because  (for  a  symmetric  well  or,  equivalently,  for  a 
two  level  system  with  zero  static  dipole  moment)  it  is  symmetry  forbidden  to  all 
orders  in  perturbation  theory  in  E0-  If  (o),Eq}  is  chosen  near  an  AD  point  the 
Fourier  transform  of  p(t)  has  no  even  harmonics  and  has  a  large  peak  at  the  low 
frequency  Q0.  The  existence  of  this  peak  is  called6'7  low  frequency  generation 
(LFG).  This  frequency  is  small,  and  can  be  continuously  timed  to  smaller  values, 
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by  changing  co  or  E0  to  get  them  closer  to  an  AD  point.  When  the  AD  point  is 
reached  becomes  zero  and  the  LFG  peak  in  p(D)  disappears  to  become  a  large 
static  dipole.  Some  of  the  work  on  this  system  assumed  that  the  electron  is 
initially  localized.  Such  a  state  is  very  difficult  to  prepare.  However,  numerical 
calculations6"7  show  that  the  same  behavior  is  obtained  if  the  system  starts  in  the 
ground  state  and  is  excited  with  a  semi-infinite  Gaussian  pulse  with  a  well 
chosen  rise  time  and  {£D,Eq}  at  or  near  an  AD  point. 

These  phenomena  were  documented  by  analysis  of  two  level  models1’2,  9,10 
and  by  numerical  calculations  with  double  wells1,2,6’7.  They  have  also  been 
observed  in  numerical  calculations  on  systems  that  mimic  an  AlAs-GaAs-AlAs 
double  quantum  well6'7;  the  analysis  of  these  results9  indicates  that,  for  the 
parameter  range  of  interest  here,  the  behaviour  of  the  quantum  double  wells 
provided  by  an  AlGaAs  superlattice  can  be  represented  rather  accurately  by  a  two 
level  model. 

In  this  article  we  investigate  how  these  phenomena  are  modified  when  the 
system  is  interacting  with  a  static  electric  field  of  intensity  Es.  The  presence  of 
this  field  introduces  a  new  parameter 

e3=  2p.12Es/h(o  (5) 

The  analog  of  the  AD  points  (  we  call  them  generalized  AD  points  (GAD))  are  the 
points  {Es,o),E0  }  for  which 

es  =  n,  (6) 

and 

e0  =  Zj(n).  (7) 

Here  n  is  an  integer  (  positive  or  negative)  and  Zj(n)  the  i-th  root  of  the  Bessel 
function  Jn.  For  parameter  values  near  or  at  a  GAD  point  the  dominant 
contribution  to  the  dipole  induced  by  the  joint  action  of  the  static  field  and  the  laser 
is: 
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(8) 


\ 

H(t)  =  p12  m,(t)  3  (es+  n)2^2  +  ^/Wn2  cosfx  Wn} 
Here  x  » t  co,  t  is  time. 


On  *  (2e/hco)Jn(e0) 

(9) 

Wn  =  [On2  +  (n-:-es)2l1'2 

(10). 

The  arguments  leading  to  equation  (8)  suggest  that  this  is  a  good  approximation 
for  p(x)  only  if  2e/hco  <  1.  The  above  equations  agree  well  with  the  induced  dipole 
calculated  numerically  only  for  sufficiently  large  values  of  e0;  the  agreement  at 
lower  values  of  e0  becomes  better  and  better  as  2e/hco  is  made  smaller. 

Note  that  in  the  absence  of  the  bias  we  have  es  =  0.  The  condition  (6)  can 
then  be  satisfied  only  if  n  =  0,  which  means  that  Wn  =  Qn  =  O0;  the  condition  v-' 
becomes  equivalent  to  Eq.  (3),  as  it  should.  When  n  is  zero,  Eq.  (8)  is  valid  for 
small  nonvanishing  values  of  es  and  generalizes  results  obtained  in  ref.  9. 

Before  we  proceed  we  comment  on  the  physical  restrictions  imposed  by 
these  conditions  and  the  systems  which  can  satisfy  them.  The  condition  es  £  1 
together  with  hco  >  2e  implies  that  the  electron  interaction  with  the  static  field 
must  be  larger  that  the  gap  2e  between  the  two  levels.  This  condition  cannot  be 
satisfied  by  atoms  or  molecules  but  are  satisfied  in  quantum  wells  because  the 
barrier  is  low  (the  electron  solid  interaction  is  weak)  and  the  system  is  large  (the 
electron  kinetic  energy  can  be  low);  perhaps  a  carefully  chosen  proton  tunneling 
system  might  also  satisfy  it.  The  condition  hco  >  2e  is  also  troublesome.  The 
system  can  be  treated  as  a  two  level  system  only  if  hco  «  eg-ej.  These  two 
conditions  require  that  e3-e2  »  £2-£i»  w^ch  is  easily  satisfied  for  a  quantum  well 
or  a  tunneling  system,  but  not  for  an  atom  or  a  molecule. 

The  results  summarized  above  are  derived  by  analysing  the  equation  of 
motion  for  the  induced  dipole.  We  show  that  p(x)  is  the  sum  of  the  result  shown  in 
Eq.  (8)  and  a  term  (2e/hco)2c(x),  where  c(x)  is  a  functional  of  p(x)).  We  then  argue 
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that,  if  the  parameters  are  sufficiently  close  to  a  GAD  point  and  2e/ho)  <  1  than 
(2E/hco)2o(t)  is  small  and  Eq.  (8)  is  a  good  approximation. 

What  follows  is  a  brief  outline  of  the  reasoning  that  lead  us  to  Eq.  (8).  We 
use  the  model  Hamiltonian: 

H  =  e{l2X2l  -1 1X1 1 )}  -  { 1 1X2 1  *  1 2X1 1  }^12[ES+E0cos(o)t)l.  (11) 

Here  1 1)  and  1 2)  are  the  ground  and  first  excited  states  of  the  electron  in  the 
symmetric  double  well  in  the  absence  of  the  laser  and  of  the  bias,  and  2e  =  e2  - 
The  zero  of  energy  is  half-way  between  the  ground  and  the  excited  state. 

The  induced  dipole  is: 

p(t)  =  P-i2  (V.t  I  ( 1 1X2 1  - 1 2)<1 1 )  I  V,t),  (12) 

where  ly,t)  satisfies  the  time  dependent  Schrodinger  equation  with  the 
Hamiltonian  (11).  By  using  standard  methods  we  derive  for  ja(t)  the  equation  of 
motion: 

dp(x)/dx=  -(2£/hco)2  2  2 

m  =  -w  k  =  -« 

Re{Jo  di^  exp[i(m+es)x]  expCiCk+e,.)-^!  pix^  (13) 

Eq.  (13)  is  used  to  search  for  points  in  the  parameter  space  where  the 
induced  dipole  may  have  unusual  behavior.  If  we  replace  in  the  rhs  p(x)  with 
we  can  perform  the  integral  and  find  that  the  terms  corresponding  to  m  =  k  =  -  es 
lead  to  secular  terms  (i.e.  they  grow  indefinitely  with  x)  which  grow  faster  with 
time  than  the  other  terms.  This  suggests  that  when  we  do  not  use  perturbation 
theory  the  terms  corresponding  to  m  =  k  =  -  es  may  be  larger  than  the  rest9. 

To  find  the  contribution  to  p(x)  from  the  terms  m  =  k  =  -  es  we  must  calculate 
them  nonperturbatively.  This  is  done  as  follows.  We  assume  that  es  -  -  n  where 
n  is  an  integer.  Then  we  write  the  rhs  of  Eq.  (13)  as  the  term  corresponding  to  k  = 
n  plus  all  the  other  terms  which  are  (2e/hco)2  c(x);  these  include  the  terms  where 
m#k  and  the  terms  where  m  =  k  but  m  *  n.  This  separation  allows  us  to  write 
Eq.  (13)  as: 
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djj-('c)/dx  =  -  ft^fj  cos[(n-+-es)(t-r1)]^('c1)}  -  (2e/tico)2c(x).  (14) 

We  can  now  use  the  Laplace  transform  method  to  “solve”  this  equation.  We 
use  the  initial  condition: 

p(t=0)  *  pl2,  (15) 

because  in  a  biased  system  in  which  p12Es  »  e  the  dipole  moment  in  the  ground 
state  (in  the  absence  of  the  laser)  is  very  close  to  p12  =  e(l  I  x  I  2).  Here  e  is  the 
electron  charge,  x  is  the  electron  position  operator  and  i  1)  and  12)  are  the 
eigenstates  of  the  unbiased  system. 

Laplace  transforming  (14),  using  (15),  solving  for  the  Laplace  transform  of 
p,  and  then  inverse-transforming  gives: 

p(x)  s  p^  pn(x)  -  (2e/hco)2  I  dx1  pn(x1)  cCx-tJ  (16) 

JO 

This  equation  is  exact,  but  it  does  not  provide  a  solution  for  our  problem  since  c(x) 
depends  of  p.  The  result  is  however  useful  because  the  last  term  at  the  right  hand 
side  is  negligible  when  (2e/ha>)  is  small;  thus  (16)  leads  to  (8).  The  argument 
suggesting  that  the  last  term  is  negligible  is  patterned  after  a  procedure  presented 
in  Ref.  9.  But  since  it  is  too  elaborate  to  be  presented  here  we  prefer  to  test  this 
suggestion  numerically. 

Eq.  (8)  predicts  that  if  eg  =»  -  n  than  p(ft)  is  dominated  by  two  modes,  a  zero 
frequency  mode  (ZFM)  and  a  low  frequency  mode  (LFM);  furthermore,  the 
equations  gives  formulae  for  the  amplitudes  of  the  modes  and  the  low  frequency 
Wn.  To  test  these  results  we  have  calculated  p(x)  numerically,  found  its  Fourier 
transform  p(ft)  and  plotted  ip(ft)l  as  a  function  of  ft.  These  plots  have  peaks 
whose  heights  give  the  amplitude  of  the  Fourier  components  of  p.  We  find  that  as 
long  as  2e/hco  «  1  and  eg  =  -  n,  the  ZFM  and  the  LFM  modes  have  the  highest 
amplitudes,  in  qualitative  agreement  with  Eq.  (8).  For  a  quantitative  test  we  plot 
in  Fig.  1  the  amplitude  of  the  ZFM  (denoted  p(0)),  the  amplitude  of  the  LFM 
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(denoted  p(Wx)  and  the  lowest  frequency  W1(  as  a  function  of  e0.  The  numerical 
calculations  were  made  for  es  =  *  1.00  and  2e/hco  =  0.129.  For  these  values  Eq.  (8) 
predicts  that  p(0)  =  0,  jifWj)  =  1  and  Wx  =  (2£yhco)J1(e0).  We  plot  these  predictions 
along  with  the  numerical  results.  The  predictions  are  reasonably  accurate  if  e0  > 
1,  and  excellent  for  e0  >  2.  The  Eq.(9),  giving  the  LFG  emission  frequency,  agTees 
extremely  well  with  the  value  obtained  numerically,  except  when  e0  is  very  small; 
Eq.  (9)  predicts  that  Wj^  0  if  e0  0,  while  the  numerical  results  lead  to  a  finite 
value. 

As  the  intensity  of  the  laser  field  is  increased  p(0)  becomes  smaller  while 
jli.(W1)  increases.  The  sum  of  these  two  components  is  close  to  the  maximum 
possible  value  jjl  12  of  the  induced  dipole.  This  means  that  the  contributions  from 
the  other  Fourier  components  are  substantially  smaller,  as  predicted  by  the 
theory.  As  e0  approaches  1  the  oscillating  part  dominates  the  dipole.  This  means 
that  practically  all  the  charge  density  is  set  in  motion  by  the  laser.  This  coherent 
motion  of  the  charge  density  leads  to  intense  LFG.  The  frequency  Wx  varies  with 
e0  (but  it  is  always  much  smaller  than  the  typical  frequencies  appearing  in  the 
Hamiltonian)  and  approaches  zero  as  e0  gets  close  to  the  first  zero  of  Jx  (which  is  a 
GAD  point).  When  Wj  =  0  the  LFM  becomes  static.  An  experiment  that  monitors 
only  the  static  part  of  p(t)  will  notice  that  the  static  dipole  disappears  as  e0  exceeds 
one  and  reappears  suddenly  when  e0  =  Zjd)  =  3.84. 

Note  that  the  analytical  results  disagree  with  the  numerical  ones  if  e0  <  1; 
for  this  values  of  e0  the  terms  neglected  in  Eqs.  (14,16)  become  important.  The 
numerical  calculations  show  that  the  values  of  e0  below  which  Eq.  iS)  is  eironeous 
become  smaller  as  2e/hc o  is  diminished. 

The  theory  makes  predictions  also  for  the  case  when  es  is  close  to  an 
integer.  In  Fig.  2  we  compare  the  theory  with  the  numerical  results  for  es=  -  0.95. 
The  agreement  is  rather  good  for  l<e0  and  excellent  for  2<e0.  As  the  laser 


intensity  is  increased  the  amplitude  ji(0)  of  the  ZFM  is  diminished  but  it  never 
becomes  zero  (as  it  does  when  es=-l).  When  e0  approaches  the  first  zero  of  Jv  p(0) 
goes  up  to  a  maximum.  At  this  point  the  electron  is  completely  localized. 

We  have  also  studied  the  case  e3=;  -  1.05.  The  equation  (8)  predicts  that  the 
results  should  be  identical  to  those  obtained  for  e3=  -  0.95,  but  we  find  this  to  be 
only  approximately  true.  Differences  between  the  calculation  for  es=  -  1.05  and 
that  for  es=  -  0.95  are  small  but  not  zero;  they  become  negligible  when  e0  >  2.5. 

The  analytical  results  are  valid  only  if  es  is  close  or  equal  to  an  integer 
Numerical  tests  show  that  they  fail  rather  badly  in  the  case  when  es=  -  1.5. 

In  the  case  of  an  unbiased  well,  studied  in  previous  work1'10,  the  ability  of 
the  laser  to  prepare  a  state  having  the  largest  possible  static  moment  (which 
means  that  the  electron  is  located  in  one  well)  is  rather  startling.  The  double  well 
studied  here  interacts  with  a  very  large  static  field  (i.e.  n12  Es>2e).  Because  of  this 
the  system  has,  in  the  absence  of  the  laser,  a  very  large  static  dipole  and  the 
electron  is  localized.  In  this  case  the  surprise  is  that,  in  spite  of  an  extremely 
strong  bias,  the  laser  manages  to  make  the  static  component  of  the  dipole  very 
small,  by  forcing  most  of  the  charge  density  to  slosh  coherently  back  and  forth 
between  the  wells. 
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Figure  Captions 

Fig.  1.  The  zero  frequency  and  the  low  frequency  components  of  the 
induced  dipole  moment  (p(0)  and  ^(W^)  and  the  lowest  frequency  Wlt  as  functions 
of  eQ.  The  continuous  lines  are  the  numerical  results.  The  values  of  WL,  p(0)  and 
pCftp  calculated  usir0  Eq.  8  are  shown  as  dots,  squares  and  x’s,  respectively.  2e  = 
12.9  cm’1,  co  =  100  cm"1,  and  es  =  1.0. 

Fig.  2.  Same  as  Fig.  1  except  that  es  =  0.95. 
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